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For high energy nucleus-nucleus colhsions, we show that 
a combination of space-momentum correlations characteris- 
tic of radial expansion together with the correlation between 
the position of a nucleon in the nucleus and its stopping, re- 
sults in a very specific rapidity dependence of directed fiow: 
a reversal of sign in the mid-rapidity region. We support our 
argument by RQMD model calculations for Au+Au collisions 
at ^/s = 200 AGeY. 
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The study of anisotropic flow in high energy nuclear 
collisions has attracted increasing attention from both ex- 
perimentalists and theorists pj-^] . The rich physics of di- 
rected and elliptical flow [p|-p^ is due to their sensitivity 
to the system evolution at early time. Anisotropic flow 
in general is also sensitive to the equation of state [p]pO| 
which governs the evolution of the system created in the 
nuclear collision. 

The collective expansion of the system created during a 
heavy-ion collision implies space-momentum correlation 
in particle distributions at freeze-out. Simplified, this 
means that particles created on the left side of the system 
move in the left direction and particles created on the 
right side move in the right direction (on average). We 
will show that the rapidity dependence of directed flow 
of nucleons and pions can address this space-momentum 
correlation experimentally. 

A sketch of a medium central symmetric heavy-ion col- 
lision is shown in Fig. |l|, from before the collision (a,b) 
to the resulting distributions of (x) and (px) shown in 
(d). In Fig. |l|a the projectile and target are shown before 
the collision in coordinate space. In Fig. |l|b the overlap 
region is magnifled and the "spectators" are not shown. 
It shows in a schematic way the number of nucleons and 
their position in the x-z plane (where x is the impact 
parameter direction). Projectile nucleons (light color) at 
negative x suffer more rapidity loss than those at positive 
X, ending up closer to mid-rapidity. Assuming a positive 
space-momentum correlation (as indicated in Fig. |l|c), 
these nucleons have negative (px), while those at positive 
X have positive (px)- This results in a wiggle structure in 
the rapidity dependencies of (x) and (px) which is shown 
schematically in Fig. 

The shape of the wiggle, both the magnitude of (px) 
and the rapidity range, depends on the strength of the 
space-momentum correlation, the initial beam-target ra- 



pidity gap and the amount of stopping. Therefore, the 
dependence of the wiggle on the collision centrality, sys- 
tem size and center of mass energy may reveal impor- 
tant information on the relation between radial flow and 
baryon stopping. In addition, it has been shown that the 
magnitude of (px) depends on the nuclear matter equa- 
tion of state llll . 
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FIG. 1. A schematic sketch of a medium central sym- 
metric heavy- ion collision in progressing time (a),(c) and the 
rapidity distribution of (px) and (x) in (d). In (b) the over- 
lap region is magnified and the "spectators" are not shown In 
these figures x is the coordinate along the impact parameter 
direction and z is the coordinate along the projectile direction 
(for a more detailed description see text). 

The above picture changes for collisions at lower ener- 
gies, where there is no clear rapidity separation between 
projectile and target nucleons at freeze-out, because nu- 
cleons cross over mid-rapidity. Moreover, when the time 
for the nuclei to pass each other becomes long relative to 
the characteristic time scale for particle production, the 
interactions between particles and spectators (so-called 
shadowing) become important. This has been pointed 
out by many people and recently in |||,|l^,|l^. This is 
consistent with the results of heavy ion collisions, in the 
2 to 158 AGeV energy range, where the experimental ob- 
served slope around mid-rapidity in directed flow shows 
a trend from a positive value at 2 AGeV to almost 
zero at 158 AGeV [|. 

Note that the change of sign of directed flow 
at mid-rapidity has been discussed for Ca-|-Ca colli- 
sions at 350 AMeV |jl^ and for Au-|-Au collisions at 
11 AGeV pT| , p^ , p^ . However, the physical origins on 
which these predictions are based are different from what 
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we discuss in this Letter. In 350 AMeV Ca+Ca colli- 
sions the wiggle originates from a combination of repul- 
sive nucleon-nucleon collisions and an attractive mean 
field. The repulsive nucleon-nucleon collisions dominate 
at the mid-rapidity region and lead to a positive slope 
for (px) versus rapidity. The attractive mean field dom- 
inates at beam and target rapidities and leads to a neg- 
ative slope. In 11 AGeV Au-I-Au collisions the wiggle is 
caused by the longitudinal hydrodynamic expansion of 
a tilted source ||l^. It has been noticed that in hydro 
calculations the wiggle only appeared if a QGP equa- 
tion of state is used. The QGP equation of state is a 
prerequisite to reach the stopping needed to create this 
tilted source. The predicted wiggle in this Letter does 
not assume a QGP equation of state. The other main 
difference is that in our prediction wc specifically assume 
in-complete stopping. 

The arguments used in this Letter which lead to a 
change of sign of directed flow at mid-rapidity are valid 
on general grounds. However, to test the picture quanti- 
tatively we study Au+Au collisions at y/s = 200 AGeV 
in an impact parameter range 6=5-10 fm, using the Rel- 
ativistic Quantum Molecular Dynamics (RQMD v2.4) 
model in cascade mode ]l9t . 

To characterize directed flow, we use the first Fourier 
coefficient P,pO|, vi, of the particle azimuthal distribu- 
tion. At a given rapidity and transverse momentum 
the coefficient is determined by vi — (cos (jj) , where (j) 
is the azimuthal angle of a particle relative to the re- 
action plane angle (x direction in Fig. |^). Similarly a 
Fourier coefficient can be determined in coordinate space, 
si = {cos(j}s) |l^, where (j)s is the azimuthal angle of a 
particle, determined from the freeze-out coordinates x 
and y, relative to the reaction plane angle. Figure |^ 
shows the RQMD calculations of vi and si for nucleons 
and pions in Au+Au collisions at RHIC energy. Indeed, 
the shape at mid-rapidity for nucleons is consistent with 
the picture described above: both vi and si show a neg- 
ative slope at mid-rapidity. For larger rapidities the si 
values leave the ordinate scale. 
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FIG. 2. RQMD calculations of vi (filled circles) and si 
(open circles) for nucleons (left panel) and pions (right panel). 

Pions are produced particles and their space-rapidity 
correlation is different from that of nucleons shown in 
Fig. Due to the asymmetry in the numbers of col- 
liding target and projectile nucleons, the pions produced 
at positive x shift toward positive rapidity. The pions 
produced at negative x shift toward negative rapidity. 
This results in a positive space-rapidity correlation with- 
out a wiggle, see Fig. ||b. (open circles). Due to the 
space-momentum correlation, the momentum distribu- 
tion tends to follow the space distribution. This leads to 
the positive slope of vi at mid-rapidity for the pions, see 
Fig. (filled circles). However, shadowing by nucleons 
is also important in the formation of pion directed flow. 
The contribution is relatively small in the mid-rapidity 
region, where in high energy nucleus-nucleus collisions 
the nucleon to pion ratio is small. At beam/target ra- 
pidity the shadowing becomes dominant, this explains 
why si has the opposite sign from vi for pions close to 
beam/target rapidity. 

In this Letter we have shown that the combination of 
space-momentum correlations characteristic of radial ex- 
pansion together with the correlation between the posi- 
tion of a nucleon in the nucleus and its stopping, results 
in a wiggle in the rapidity dependence of directed flow 
in high energy nucleus- nucleus collisions. Moreover, the 
amount of stopping and the space-momentum correla- 
tion depend on the equation of state and this affects the 
strength of the wiggle around mid-rapidity . Finally, 
because the wiggle appears at mid-rapidity, it is accessi- 
ble by the current SPS experiment NA49 and the near 
future RHIC experiments. The study of its dependence 
on collision centrality, system size and the center of mass 
energy may reveal important information on the relation 
between collective radial flow and baryon stopping. 
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